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WHY DO WE NEED 
PROPULSION TESTING* 


AND 
BOOKKEEPING 


Presented by 
Bob Berrier 
4/17/17 
* Care should be exercised as this nomenclature has been used 
for Inlet tests, nozzle tests, engine tests, and PAI tests 


AFT-END DRAG AT SUBSONIC SPEEDS 


Later 1960s 1970s 1980s 


Early 1960s 


INLET/AIRFRAME INTERFERENCE EFFECTS 


Le Interference drag 


Design mass flow 


Cp ac 


Inlet mass flow 


TYPICAL FIGHTER PROPULSION INSTALLATIONS 


TYPICAL BOMBER PROPULSION INSTALLATIONS 
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TYPICAL TRANSPORT PROPULSION INSTALLATIONS 


PROPULSION INTEGRATION TASKS 


Basic Research/Parametric Design Data 
o Nonaxisymmetric/Multifunction Nozzles 
o Tail interference 

o Theory/CFD validation 


Product Improvement 
o F-111, F-14, F-15, F-18, B-1B, B-2 


Evaluation of New Concepts 
o Distributed Propulsion 
Blended Wing Concepts 
Boundary Layer Ingestion 
Low Boom Supersonic Config. 


ORS ae ©. 


Flight Performance Evaluation (Valid bookkeeping system that 


is agreed to by everyone is essential) 
o Corrections to Flight Aerodynamic and Engine Data Packages 


TYPICAL BOOKKEEPING SYSTEM AND REQUIRED MODELS 
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SUCCESSFUL BOOKKEEPING RESULT — THRUST/DRAG CLOSURE 


TOTAL AIRCRAFT SYSTEM PERFORMANCE SIMULATION 


Historical Updated Updated Model ta 
Drag Drag Installation Full Scale Final 
Polars Polars Drag Corrections Corrections 
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AIRFRAME Pre-Fight]  /Ciigh Compare A/F 
PROPULSION — Friction = a rel Test a a |.) Prapu lsion 
: : % Installed Full-span Drag a 
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Propulsion 


Historic leclated Instrumentation 


Nozzle Suppression Installed Corrections 
Data Suppression Eic. 
Data 


solated Inlet Drag 


Nozzle Tests Tests or Analysis Model to Data Final 
* Suppression TPS Inet Mir. Full Scale Prvestine 
- Nozzle Coeff. v5 Corrections 
Flight Inlet Mfr. 


CNGINE PERFORMANCE SIMULATION 


ATF Perf. incl. Update ATF 
eee Flight Suppression, Performance Incl. 
Scat C,, Update Installed Suppression 


Test of Flight 
Test Engine 


TYPICAL PROPULSION 
TESTS/MODELS/RESULTS 


WHAT ARE THEY USED FOR? 


PROPULSION COMPONENT TEST CONTRIBUTIONS 


INLET SYSTEM ENGINE 
¢ Part of AC bookkeeping ¢ Where full-scale thrust 
system | value comes from 
* Provides | ¢ Separate bookkeeping 
— Pressure recovery | system 
— Distortion ¢ Uses info from inlet & 


exhaust tests and vice- 
| versa 


— Operating geometry A 
— Spillage drag A 


FOR AC PERF, ALL | 
PARTIES MUST BE IN 


AGREEMENT | 


EXHAUST SYSTEM 
¢ Part of AC bookkeeping 
| system 
| * Does not provide full-scale 
thrust (not scalable) 
¢ Provides 
| — C, (Nozzle efficiency) 
| - C, (Discharge coeff.) 

— Sting and Distortion A 
| — Operating geometry A 
— Operating NPR A 
— Engine suppression 
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ENGINE MODEL 


(caculated)  & TYPICAL ENGINE BOOKKEEPING 
SYSTEM - ONE OF MANY THRUST 
METHODS FROM SAE AIR 1703 
AND SAE AIR 5450 
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DRAG COMPARISON AT M=0.85 (7) 
Ry/FT = 3.85 x 108 


O 16FT. (FEB.1971) : 
O 8FT.(AUG.1969) Drag Polar 


¢ Fixed variables 
— Mach No. 
— Inlet Capture Ratio, Aj/A; 
— NPR (Flow-Thru) 
— Roll/Sideslip Angles 
— Control Settings 
— Nozzle/Inlet Ext. Geometry 
-Ry 
— Weapon Load 


Lift Coefficient, Cy 


TYPICAL DRAG POLOR CORRECTION ON F-111 


Over Wing Fairing 


yn A =26° 
Flexibility Base M=0.86 
Secondary Air Spillage q = 450 lb/ft? 
Proturberence 

Roughness 


Reynolds Number 


Laminar to Turbulent 


iz 
Themed Full Scale 
Wind Tunnel ———— Estimate 
.008 .012 016 .020 024 .028 .032 


Drag Coefficient, Cp 


FIGURE 72 - ADJUSTMENTS MADE TO F-111 WIND-TUNNEL DRAG (AFFDL-TR-78-100) 
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SPILLAGE DRAG CORRECTION 


(=) _ Captured streamtube freestream area 


indet projected capture area The throttle-dependent inlet spillage drag 
force increment is defined as the change 
in aircraft drag force resulting from the 
difference between operating and 
reference inlet mass flow ratios 


Ac 


S ref = constant 


Inserts 


Incremental inlet spillage drag coefficient, 


0 : 1.0 
Mass flow ratio, Ag/Ac 
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EXTERNAL INLET PERFORMANCE CORRECTIONS 


F-15 INLET SYSTEM 


hh 
Second and Third Ramp Bleed Exits 
Sip Za 
Bypass Door (Variable) <4 
Slot Bleed/Bypass 
Exit (Fixed) 
Diffuser Ramp 
—™ = — 
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me q 
Second Ramp 
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inlet aerodynamic force 
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Moment 
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Inlet capture ratio, Ao/Ac 
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INLET MODEL 


Inlet performance 


Recovery 


Distortion 
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Inlet capture ratio, Ao/Ac 
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INTERNAL INLET 
PERFORMANCE 


EXPERIMENTAL COMPUTATIONAL 


M = 0.25 
Re,» = 6.9 x 10° 
W2C/A, = 294 


io 
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FAN FACE DISTORTION MAPS 


ms 


M = 0.83 
Re,» = 13.9 x 10° 
W2C/A, = 29:8 

a54 


M = 0.25 
Re, = 6.9 x 10° 
W2C/A, = 32.0 
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M = 0.83 
Re,» = 13.9 x 106 
W2C/A, = 36.8 
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EXHAUST SIMULATION 
MODEL 


F-15 EXHAUST SIMULATION MODEL 


Thrust coefficient, Cy 


INTERNAL NOZZLE PERFORMANCE 


Ag/Ag 


Nozzle pressure ratio, Pyp/P.o, 


Discharge coefficient, Cp 


Nozzle pressure ratio, Pyp/Pe, 
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Change in nozzle 
discharge coefficient, ACy 


NOZZLE FLOW COEFFICIENT SUPPRESSION 
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0.8=M \— Descent 


tet 1.4 1.6 1.8 2.0 2.2 2.4 
Nozzle pressure ratio, Py47/Po, 
O90916 


STING AND DISTORTION CORRECTION MODELS 


STING AND DISTORTION CORRECTIONS 


Sting and Distortion Increments, Cruise Nozzle, Swept Strut 4 


Flow-through pressure ratio 


a= (0° b, = (0° 
. 0060 : LRC propulsion model 
0040 MCAIR propulsion model 
0020 
ACD 

0 
Performance baseline 

—*. 0020 
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JET INTERFERENCE CORRECTIONS 
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JET INTERFERENCE INCREMENTS: SWEPT STRUT 


Propulsion Model 
= = Oo 
a= 0° O, 0 


ae Cy engine PR ~ “D, flow-through PR 
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CRUISE NOZZLE 


O LRC propulsion model 
O MCAIR propulsion model 


LOW MACH A/B NOZZLE 


SOME REFLECTIONS ON 
AERODYNAMIC FORCE & MOMENT 
TESTING 


Francis J. Capone 
NASA Retired 


April 17, 2017 


AERODYNAMIC FORCE & MOMENT 
TESTING 


OUTLINE 


MODEL SIZING/INSTALLATION 
BALANCES 

— Installation 

— Calibration 
STING DESIGN 

— Strength Characteristics 

— Sting Divergence 

— Interference Effects 
WEIGHT AND ATTITUDE TARES 
MODEL ATTITUDE 

— Importance 

— Basic equations 

— Tunnel Flow Angularity 
OTHER 
CORRECTIONS/ADJUSTMENTS 

— Base Pressure 

— Nacelle/Internal Duct Drag 
BOUNDARY LAYER TRANSITION 
OTHER CONCERNS 


MODEL SIZING/INSTALLATION 


HIGH-SPEED MODEL SIZING CRITERIA 


¢ Model Cross-Sectional Size Dictated By Blockage (A,,o¢e/Atunnel): 
— Subsonic speeds- M<1 
— Supersonic speeds - 1.00<Ms 1.04 


¢ Model Length Dependent On: 
— “Calibrated” test-section length at subsonic/supersonic speeds 


— Boundary reflected disturbances at supersonic speeds 


¢ Model Volume: 


— Dependent if test section Mach number gradient sufficient to require 
buoyancy correction to drag coefficient (accepted gradient dM/dL < 0.0006) 


— Buoyancy correction = f(Volume)(dp,,/dL) 


MODEL SPEED RANGE CONSIDERATIONS 


¢ Model Flow Field/Wall Interference Effects Considered for the 
Following Mach Number Ranges: 
— M< 1.00 - Subsonic/Transonic blockage 
— M= 1.00 to 1.05 - Supersonic blockage 
— M> 1.05 - Shock reflections 


¢ Results Presented Pertain to Wind Tunnels With Slotted Walls 


¢ Blockage And Shock Reflection Criteria Derived From Specific 
Studies That Investigated These Problems 


SUBSONIC BLOCKAGE INTERFERENCE 


Tunnel Blockage 
»—= = 16 FIT 0.068% 
"== 8FIT 0.280% 


0.94 0.96 0.98 1.00 1.02 
Mach number 


¢ Shown is a comparison of experimental to flight drag rise characteristics for the same 
supercritical body of revolution. 
Because of the very poor correlation of these results, a detailed study of blockage effects of 
thirteen bodies of revolution with four different profile shapes was initiated in both the 16-Ft 
Transonic Tunnel (16 FTT) and the 8-Ft Transonic Pressure Tunnel (8 FTT). 


All four body types exhibited the same drag characteristics as shown above, i.e., drag of the 
body was significantly reduced as blockage ratio was increased. 


Couch and Brooks: Effect of Blockage Ratio on Drag and Pressure Distributions for Bodies of Revolution 
at Transonic Speeds. NASA TN-D7331, 1973. 


SUBSONIC BLOCKAGE EFFECTS ON BODY 
PRESSURE DISTRIBUTIONS 


a ~.. Blockage pressurizes 
| Lo — over 40% of the body 
-0.0 


i gee = 


== 8FTT 0.140% 


¢ One conclusion from NASA TN D-7331 was 


there was only one effect of wall interference on the model surface-pressure 
distributions obtained for a given model shape at different values of blockage ratio. An 


increase in blockage ratio for Mach numbers greater than 0.96 caused a positive 
increment of pressure to occur on the model” 


¢ The net result of this increase in pressure was a decrease in drag. 


SUBSONIC BLOCKAGE CRITERIA 


Experimental Data, TN D-7331 
2.0% Model, 0.156% Blockage 
3.0% Model, 0.350% Blockage 
3.5% Model, 0.480% Blockage 
4.0% Model, 0.630% Blockage 
5.0% Model, 0.980% Blockage 


0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 
Blockage ratio 


* NASA TN D-7331 also concluded that: 


— “increasing the blockage ratio above approximately 0.0003 produced a premature, positive deviation, or 
transonic creep, of the drag curve from the trend of the subsonic data. Since the Mach number at the 
initiation of the transonic creep agrees with the calculated Mach number for choked flow in a solid-wall 
tunnel, transonic creep may be the first indication of significant wind-tunnel wall interference near a 


Mach number of 1.0.” 


¢ This was the first investigation to indicate that the flow relief afforded by the test section slots 
- designed according to previously accepted criteria for interference-free subsonic walls - did 
not appear to be sufficient to avoid significant interference of the walls with the model flow 
field for Mach numbers close to 1.0. 


WALL INTERFERENCE IN SLOTTED WIND TUNNEL AT 
MACH NUMBERS FROM 1.00 TO 1.04 


al 
At Vy \ V \ 


i eo p/p, 


“0 1 2 3 4 5 6 7 8 9 10 11 12 


¢ Wall interference effects at these Mach numbers are dominated by reflection of the cone- 
shoulder expansion field being reflected back to the model as compression waves. 


The location of the reflected compression waves is dependent on the model blockage 
ratio. 


¢ Capone and Coates: Determination of Boundary-Reflected-Disturbance Lengths in the Langley 16- 
Foot Transonic Tunnel. NASA TN D-4153, 1967. 
¢ Similar results reported in NACA TN 4233 and AEDC-TR-59-12 (perforated walls). 


SHOCK REFLECTION CHARACTERISTICS FOR 


SLOTTED TUNNELS 
— m —— | ae po 
a | oA Ea 
iN , fatto p th 
a ais 
fas SPT Dea 
ar 

pa | Reflected Disturbances 

0 1 2 3 4 5 


8 9 10 11 12 
Wall interference effects at Mach > 1.05 are dominated by impingement of 


boundary-reflected-disturbances (reflected bow shock) on the model 


Shock reflection length, in 


SHOCK REFECTION LENGTH CRITERIA FOR 
SLOTTED TUNNELS 


Measured shock refection lengths for 16 FTT Estimated shock refection lengths for NTF 


= 
= 
Dd 
< 
2 
< 
= 
eae 7 
Interference free 2 
= = NACA 1389 9 
.04 1.06 1.08 1.10 1.12 1.14 1.16 1.18 1.20 1.04 1.06 1.08 1.10 1.12 1.14 1.16 1.18 1.20 
Mach number Mach number 


Results indicated that the stronger the bow shock wave, the shorter the maximum interference-free 
model length. 


Shock reflection lengths determined from the Mach angle are much longer than those measured. 
The dashed line represents the average shock reflection lengths measured in the original 12-sided 8 


Ft. tunnel (NACA Rpt. 1389) using a body with a nose half angle approximately 10°. 


Results for additional configurations can be found in NASA TN D4152 and NASA/TM-1998-208723. 
The body used in the latter reference approximated the fuselage of a supersonic transport. 
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Shock Reflection Probe (SRP) Shock Reflection 
L hs For 16 FTT & BTWT 
none rere Shock strength for BSWT 


COMPARISON OF REFLECTED SHOCK CHARACTERISTICS 
FOR SLOTTED AND PERFORATED TUNNELS 


el nea eae iris “Reflected shock-16FTT 0.05 
0.054 | 0.04 
0.0 
CP 0.001 : aC, 
! 02 
aalt —— 16 FTT - Octagonal o0 
-0.0 BTWT-Reflected shock/ — — BTWT- Rectangular 0.01 
sidewalls , 
-0.1 0 | T 0.00¢4A?4-+-_+—__+— 
0 10 20 30 40 50 60 70 0 1 2 3 4 5 6 
Body length, in. % Porosity 


A body of revolution with a nose shape that produces a bow shock of approximately the same 
strength as a realistic HSCT configuration was tested in the: 

— Langley 16-FTT slotted octagonal test section 

— Boeing BIWT ~ 8x12 slotted rectangular test section 

— Boeing BSWT 4 foot perforated test section insert 
A conclusion from AIAA-94-1935 was “ The long standing rule to not test when reflected 


shocks impinge on a model is valid and appropriate for slotted-wall test sections; slotted 
walls do not provide significant attenuation of shocks.” 

Another conclusion from AIAA-94-1935 was “Any force, moment, and pressure data 
obtained at low supersonic Mach numbers should be viewed with caution when shock 
reflection occur on a model, since shock reflections significantly affect the flow.” 
REFERENCES: NASA/TM-1998-208723, October, 1998. AIAA-94-1935, June, 1994. 
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BLOCKAGE & SHOCK REFLECTION EFFECTS ON 
SHUTTLE AXIAL FORCE 


Blockage Ratio Model Scale Tunnel MaxM.. 
O 0.0008 2 % LaRC 16 Ft 0.95 


L) 0.0053 5 %o LaRC 16 Ft ~0).90 
~ 0.0114 5 % ARC 11Ft 2? 


Max M based on NASA TN D-7331 


Blockage effects 


16 
12 
C 
A 08 Reflected bow 
shocks 
04 


Range of blockage/shock 
reflected data 


© 6 7 8&8 99 10 1.1 #12 #13 «14 


BLOCKAGE & SHOCK REFLECTION EFFECTS ON 


SHUTTLE PITCHING MOMENT 


Blockage Ratio Model Scale Tunnel Max Minter free 
© 0.0008 2% LaRC 16 Ft 0.95 


L] 0.0053 5 % LaRC 16 Ft ~0.90 
© 0.0114 5 % ARC 11 Ft ae 


Range of blockage/shock 
reflected data 
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TEST SECTION “LENGTH” 


¢ Once a model length is known, it is necessary to take a look at 
the length of the test section available at the highest Mach 
number planned for the test. Is it adequate for the model to be 
positioned properly longitudinally? 


¢ For example, test section length was dependent on Mach 
number in the 16-Foot Transonic Tunnel 
— Maximum model length = 25 feet for M < 1 
— Maximum model length = 8 feet for M = 1.20 
— Maximum model length = 10 feet for M = 1.30 


¢ If the model is outside calibrated test section, are buoyancy 
corrections required? 


FORCE BALANCES 


SOME GENERALITIES ABOUT FORCE BALANCES 


Langley 
— Fabricated from a single piece of metal 


Flexible for greater output resolution 

Temperature compensated for modulus change 

Combined load calibrations 

Performance documented over time by instrumentation technicians 


Task 
— Assembled from parts (floating frame construction) 
More rigid than langley balance for same load range 
Often not compensated for modulus change with temperature 
Some contain thermocouples for computing modulus temperature effects 
Can have internal mechanical stops (for safety) unknown to user 


Should be completely calibrated at Langley over planned test range unless 
customer supplied documentation satisfactory 


¢ Other 
— Boeing balances usually perform well 
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TYPICAL LANGLEY FORCE BALANCE 


QS 2.000 Dia.t 2.060 Dia! Keyseat 0.2500 x 0.250 deep x 1.250 long 
1.4160 max dia. 
0.3750 dia. thru hole sing Wx per ieck Gh ac Nabe 


0.6250 dia. dow! nole 
1.562- 16 thread 


1621 


rear a 
[Normal [S000 T-Nomal [S000 


[Axial | 5001S [Axial | 8001 
[Pitch [70,000 1n-16 [Pitch [70,000 in-16 
[Roll _|_7500in-10_ [Roll __|_7500in-1_ 
[Side [18001 


PREPARATION FOR TESTING - 
BALANCE INSTALLATION 


¢ Balance Fit: 
— Check fit of balance to model balance-block 
— Provide custom made tight fitting balance-to-model block pin(s) 
— Check taper surface contact between balance and sting tapers 
— Provide tight fitting roll key for use in balance-to-sting tapers 


— Check to see that the balance-to-sting roll key length or height does not prevent balance taper 
from seating 


— Use AMS units to document quality of fits of assembled support hardware 
— After all tubes and leads passed through sting, plug sting passage with clay 


¢ Pre-Model Continuity Checks: 
— Before model buildup when balance taper and mini-plug are accessible, hand load all balance 
components through the DAS 


¢ Balance/Model Loadings In Tunnel/Model Buildup Bay 
— Determine balance/sting deflection constants (backup/alternate to model attitude transmitter 
measurements 


— Determine restraint effects on balance readings due to tubes, leads, etc. 
— Assure adequate metric/nonmetric clearances under load (foul-free setup) 
— In-tunnel check loads to verify force measuring system as installed 
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ONE SOURCE OF ERROR IN a MEASUREMENT 


Roll Angle Due to Slop in Taper Key 


Taper diameter (D) 
\ | K— peo—> 4 


For a given key slop of 0.001" 
this plot shows the angular slop 
as a function of taper diameter 


Roll slop, deg 


Taper diameter (D) at keyway, in. 
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BALANCE FORCE CALIBRATION 


Double knife edge 


Stin 
xX g 


Bubble level 


Balance fixture Balance pitch center 


Calibration rig 


Weights 


NOTE: 
1. Level balance fixture between loadings 


2. Roll sting and balance 180° at sting coupling 
to get negative force calibration 
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PREFERRED BALANCE MOMENT CALIBRATION 


Moment = 1(Aweight) 


-— 1 
Bubble level 


Calibration rig 


Double knife edge 


Zo 


Double knife edge 


Sting 


Balance pitch center 


Balance fixture 


NOTE: 


1. Moment obtained by shifting weights 


2. Level balance fixture between loadings 
3. Roll sting and balance 90° at sting coupling 
to make yaw calibrations 
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STING DESIGN 


STING DESIGN - STRENGTH 


Max stress point 


4340 Steel heat treated 
to 1&0 KSI ultimate 


Allowable Stress 


p06 | CUO 19,634 45,000 
a ae 
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STING DIVERGENCE CRITERIA 


Fry AFy 
Slope = Wa = Ena (0)S #/° 


Fry aft 2 \ 
ae) (0.4)AP 


A 
Divergence criteria < 
g Aa 


P AP 
Slope = AO #/ 
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STING DIVERGENCE - ANOTHER METHOD 


¢ Sting will diverge when 


— d/L = dwido = 2El,,/ 12 


— where 
1 length 
L Lift 


A angle of attack 
E modulus of material 
Ix moment of inertia at critical section 


— now: A/L = Cha (Q/Pt) Steg Pt 


— where 
Cie lift curve slope in radians 
q/p, ratio of dynamic pressure to stagnation pressure 
Srer Pt is the product of reference area and stagnation 
pressure 
— Therefore: For sting divergence 


Dt e (2El,, / 12){1/[CL a (q/p;) Stet p]} 


— For safe practice it is recommended that p, for sting 
divergence be at least twice the test p, 


Critical section 
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STING INTERFERENCE 


Misys 


Sting flare 


Factors Affecting Sting Interference 


— d/D,,5- Ratio of sting diameter to base diameter 


- Dbase Length of cylindrical portion of sting 
- 6 Sting flare (taper) angle 


Sting Interference “Rules of Thumb” 
— With @ = 8°, 1/D,4¢ should be = 4 if M21 
— For large 6, 1/D,,;- Must increase 


— For d/Dyase £ 0-9 there are only small effects on 
“forebody (boattail) drag” (large effects on base 
drag) Cor = Co total 7 Cp base 


OTHER STING DESIGN CONSIDERATIONS 
STING LOCATION EFFECTS ON BASE PRESSURE COEFFICIENT (M = 0.9) 


D | Ds 
X 
—<l eee  OT. 
Taper angle 
0.10 
sf Config. D, in. 
0.08 \ O 15° taper 1.5 
' O 25° taper 1.5 
IN A 15° taper 1.0 
0.06 AC,» = Coy (x/Dg) - Cp p(x/D, = 6.9) 
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OTHER STING DESIGN CONSIDERATIONS 


REFLECTED SHOCK LOCATION EFFECTS ON BASE PRESSURE 


Y 


s3252 222 pad 
77 EBB: 
2F = cee 
{ 
+ $ 


«(Cc .) 


p,b interference ~ Cy) 


p. dD no interference 
M#1.53, R= 35x 10°, ref. 14 
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OTHER STING DESIGN CONSIDERATIONS 
STING-TO-BASE DIAMETER EFECTS ON AFTERBODY DRAG COEFFICIENT 


0 
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


b 
Variation of jet-off afterbody drag with sting size 


Reference: AEDC TR-80-8 


Cp based on body cross-sectional area 
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OTHER STING DESIGN CONSIDERATIONS 


MODEL BOATTAIL C,,, AND CRITICAL STING LENGTH RATIOS 


-d/D, = 0.87, 6 = 12° 


M140 


o Yp,+0 
a Yo, 5. 
© Vo,+ 85 
& 
a 
ape 
oe 
i?) 


8 


go? 
690° Ctosesesaganesee® 


“a 
o, 
on 


(Yop), for Cp, 


o 

L=) 

. te ed 

2 ee 
SERRERS. = 2B 


~— 
PT TELL Ce 


Reference: RM A57109 
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OTHER STING DESIGN CONSIDERATIONS 
MODEL BASE PRESSURE COEFFICIENTS FOR VARIOUS STING CONFIGURATIONS 


° 


D = 1.0-in. Model Diameter 
Ds = Maximum Sting Diameter, 1.5 in. 


-0.24 


a x,/D  Xo/D  x4/D,  ¢,, deg 5, deg 
nial, ro) 8.4 10.3 6.9 0 15 
“0.16; — ra) 0.5 4.0 2.7 2.0 15 
Cp, 4 2.0 5.5 3.3 2.0 15 
-0.12 © 0.5 5.5 3.7 2.0 25 
7) 2.0 25 


Reference: AEDC TR-80-8 
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A REFRESHER ON WHAT YOU SHOULD 
ALREADY KNOW 


NF uncor Raw —T Readings 
10 Apply Balance Interactions 
Using Initial Tare Loads 
AF Uncor I 


Apply Attitude Tares 


Body Axis Components Determine a and B 
Rotation & Translation <=] (Apply tunnel angularity 
Balance to Model Axis corrections) 


V 


Apply corrections for: 
Base Pressure 
Nacelle Internal Forces & Moments 
Buoyancy 


Determine Coefficients About: 
Body Axis 
Stability Axis 
Wind Axis 
Alternate Axis 
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INITIAL MODEL WEIGHT 
AND ATTITUDE TARES 


DEFINITIONS 


¢ INITIAL WEIGHT TARES: Forces and moments exerted on 
balance due to weight of the model when balance is is at attitude 
for recording wind-off data zeros 
— Used only in calculation of balance interactions 


¢ ATTITUDE TARES: The changes in forces and moments felt by 
the balance due to model weight as model (and balance) attitude 
is changed relative to the wind-off attitude 
— Each data point must be corrected for attitude tares to obtain the 
appropriate aerodynamic data 
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MODEL WEIGHT RELATED IN SIMPLIFIED TERMS 


Initial Weight Tares Attitude Tares 
An interaction matrix relative to the 
unloaded (bare) balance condition 
Is obtained from the balance 
calibration. Then, in use, the 
balance is treated as a differential 
force and moment transducer. 


NFO(I) = -W AFTARE = W(sin 6- sin8,) 
AFO(I) =0 NFTARE = -W(cos 6 - cos8,) 
PMO(l) = -Wx PMTARE = AFTARE(z) + NFTARE(x) 


DETERMINATION OF MODEL INITIAL WEIGHT AND 
ATTITUDE TARES 


To Calculate the Initial Weight Tares and Attitude Tares the Following 
Must Be Determined: 


1. The weight (W) of the model. 


2. The location of the model reference center of gravity relative to the 
balance moment reference. (x, y, Z). 


This Can Be Done By: 


1. Use of scales, etc to weigh the model and to determine the center of 
gravity location prior to the wind tunnel test. 


2. Use the balance and attitude variation in the tunnel to “weigh” the 
model (the easiest and most common approach) 
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IN-TUNNEL DETERMINATION OF MODEL WEIGHT AND 
ATTITUDE TARES 


¢ METHOD 1: Pitch model through maximum possible of support mechanism angle 
recording data at discrete angles. The data system then recognizes a run code and 
automatically uses AF output to calculate model weight. Appropriate moment outputs 
(depending on model roll) are used to determine moment tares and x, y, and Z. 


¢ METHOD 2: Roll model through 270° using strut roll mechanism recording data at each 
90° roll increment. The data system then uses NF and SF output to determine weight. 
PM or YM output are used to determine pitch (or yaw). RM output is used to determine 
roll tare. 


¢ The results of either method should be examined critically before proceeding to the test 
runs since unique circumstances exist that can produce bizarre results (c.g. 24 inches 
from the model base) when balance resolution is too small or the distance between the 
balance moment reference center and the model c.g. are very small. Results of either 
method or the two methods combined provide a complete set of initial weight and 
attitude tares that are automatically input for subsequent tunnel runs. 


MODEL ATTITUDE 


ATTITUDE DETERMINATION NOT SO SIMPLE 


“, 
“, 
[/ “~ Op >: 
a 0 Q 


Qlup 
ae Wind 


8, Strut or support system pitch angle 

0. Euler pitch angle to account for knuckle or offset sting 

0) Wind off attitude of balance 

8, Euler pitch rotation between balance and model reference axis 
8 Model Euler pitch angle 

Qyp Upflow angle 

a Model angle of attack 


F. J. Capone, 03/13/2008 a 


DETERMINE MODEL ATTITUDE IN TUNNEL 
THEORY AND BACKGROUND 


To reduce wind tunnel aerodynamic data to 
coefficient form, it is necessary to determine x, y, Z 
components of various velocity vectors, force vectors, 
and moment vectors in any of the following axis 
systems: 


Wind Axis 

Stability Axis (not shown) 
Gravity Axis 

Balance Axis 

Body Axis 


Missile Axis (not shown) 
42 


DETERMINE MODEL ATTITUDE IN TUNNEL 
USE EULER ANGLE THEORY 


¢ Use of Euler angle theory to rotate to axis systems 


¢ The components of a vector in one axis system are rotated to 
other axes systems for some data reduction systems. 


¢ That is 
— Any axis system can be aligned with any other axis system by a 
translation and three rotations, providing that the rotations are a yaw, 
pitch, and roll. 


— Any series of consecutive rotations are equivalent to some yaw, 
pitch, roll rotation. 


— The resulting rotation matrix is orthogonal 
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DETERMINE MODEL ATTITUDE IN TUNNEL 
AXES DIFINITIONS 


Wind Axis 
Zw 
sh Gravity Axis 
‘3 ei eo 
ae 
- a Yg we i Approximate direction of V.. 
V Xo 
a 
a 
Body Axis | Balance Axis 
Weight Zon 
Ze a 
Ybal NF 
x 
= Xp SF (ae 


Balance Pitch 
Center 


ROTATION OF FORCES 
GRAVITY AXIS SYSTEM TO BALANCE AXIS SYSTEM 


First Rotation - Yaw 


X' = X,cosy — Y ,siny 


ay . Y' = X,siny + Y,cosw 
g Lia Oe 
Or 
y' X'| |cosy -siny 0 ||X, 
Y'|=|siny cosy 0 |/Y, 
x 7 0 0 IZ 
X or 


F’=T,F, 


F. J. Capone, 03/13/2008 = 


ROTATION OF FORCES - Cont. 
GRAVITY AXIS SYSTEM TO BALANCE AXIS SYSTEM 


Second Rotation - Pitch 


Z' X”" = X'cos8 — Z’sinO 
Vhs y"=y' 
Z" = X'sin8 + Z'cosO 
or 
X"| |tcosO OQ —sinO | |X’ 
ie ba -0 1 0 |ly 
ZL sn8 0 cos0 Zz 
x" or 


F" = T oF’ 


F. J. Capone, 03/13/2008 46 


ROTATION OF FORCES - Cont. 
GRAVITY AXIS SYSTEM TO BALANCE AXIS SYSTEM 


Third Rotation - Roll 


oi Xpal — pe 
Lyal Ypat = Y'coso — Z” sino 
Za = Y'sind + Z"coso 


Or 


: Xpal l 0 0 “i 
Y Ypar f= |O9  cosm —sino 7 
Lal QO sind cosd . 

Y pal os 


Meats Tyke 


F. J. Capone, 03/13/2008 47 


COMBINED ROTATION OF FORCES 
GRAVITY AXIS SYSTEM TO BALANCE AXIS SYSTEM 


Foal = Tie 
Ea Ty (TF ‘) 
oh Ty ToT F,) 


Therefore: 


F, 1 0 QO |{cos0 0 -sinO 
F,| =| 0 coso -sing || 0 | 0 


cosy -siny 0 
siny cosy 0 


Fo }oa | sind coso |binB O- cosO })0 0 


(Roll) (Pitch) 


F. J. Capone, 03/13/2008 


(Yaw) 


| 


F 


Xx 


Fy 


F 


Z 
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COMBINED ROTATION OF FORCES 
GRAVITY AXIX SYSTEM TO BALANCE AXIS SYSTEM - Cont. 


Now: 


Substituting and performing the matrix multiplication gives: 


Py cosOcosy —sinycos@ -sinO 0 
Fy = |-sindsinOcosyw+cososiny —§ sindsinOsiny+cosdcosy — -sindcosO 0 
Fo Iba | COSPsinOcosytsingdsiny —cosdsinOsiny+sindcosy — cosdcosO —W 
Therefore: Fy bal = WsinO 
Fy bal = Wsingcos® 


F, bal — ~>WeosocosO 


F. J. Capone, 03/13/2008 49 


CALCULATION OF MOMENTS PRODUCED BY MODEL 


WEIGHT VECTOR IN BALANCE AXIS SYSTEM 


Model TT) 


ref center Mpal 


Ibal 


F. J. Capone, 03/13/2008 


Pitch: 
M,,) — XF, + ZF, 
Myal = -WxcosocosO + Wzsin® 


Yaw: 
Dpal xP, a yo, 
Nyal = WxsindcosO + WysinO 


Roll: 
Lal = -yF, Bs ZF, 
lhal = Wycosocos8 + WzsindcosO 
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TARE COMPONENTS OF THE MODEL WEIGHT VECTOR IN 
THE BALANCE AXIS SYSTEM 


To Summarize 


The tare components of the model weight vector in the balance axis system 
are: 


AF =F, 5a, = Wsino 

SF =Fy pa, = Wsingcosé 

NF =F,,,, = -Wcosocosé 

PM =m, 1 = -Wxcosdcosé + Wzsind 
YM =n,,, = Wxsinocosé + Wysiné 


RM =l,41 = Wycosdcosé + Wzsindcosé 


Now: The initial weight tares for the model/balance are computed using the 
above equations and the values of 8 and 9 at data zero (0) and 4) 


F. J. Capone, 03/13/2008 ot 


CALCULATION OF ATTITUDE TARES 


AP apa Fyal.0 


Therefore: 


ASF =F),4. = W(sindcos8—-singycos8y) 

ANF = F,441 = —W(cosdcos0—cosogcos8y) 

APM = my, = —Wx(cosdcosO—cosdycos0y) + Wz(sinO—sin8,) 

AYM = ny, = Wx(sindcos0—singycos0,9) + Wy(sin8—sin0) 

ARM = _l,,, = Wy(cosdcos0—cosdycos8y) + Wz(sindcosO—singycos8,) 
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DETERMINE MODEL ATTITUDE IN TUNNEL 


AXES ROTATIONS 


¢ Rotation From Wind to Gravity Axes 
— Following order will always be used: 


l. 
2; 


Wu 
0, 


sideflow 
upflow 


¢ Rotation From Gravity to Balance Axes 


— Following order should be used** (Note: Start at support strut and work forward) 


l. 


ig aa el 


strut pitch angle 
strut roll angle 
yaw knuckle 
pitch knuckle 
pitch deflection 
yaw deflection Order not important since these angles are small” 


roll deflection 


Order dependent on installation 


*For sweep tests, 8, should be first, for b sweep tests W, should be first, however 


**A notable exception to the order given is when “angle of attack” indicators in model are used 
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DETERMINE MODEL ATTITUDE IN TUNNEL 
AXES ROTATIONS - Cont 


¢ Rotation From Balance to Body Axes 


— Following angles will be used in order necessary: 
¢ 6, pitch 
° pe roll 
° Ug yaw 


¢ Rotation From Gravity to Balance Is Required for 
Weight Tare Calculations 


¢ Rotation From Wind to Body Is Required to Determine 
a and B 
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DETERMINE MODEL ATTITUDE IN TUNNEL 
CALCULATION OF Q AND B 


o = tan-'(w/u) and B = sin'(v/V..) 


Where u, v, w are x, y, Z components of V., in body axis system. 
Therefore, a rotation from wind axis system to body axis system is 
required. 
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FLOW ANGULARITY DETERMINATION 


. Starts long before the test by making sure the model can be leveled precisely upright and inverted 
whenever necessary during the test. (Provision should be made for leveling when the model is at +90° 
for tests requiring sideslip testing). Make leveling fixture identical weight to any model component 
removed to accept fixture. 


. Level model inverted and run a baseline configuration at 8= -2° to 2° at 0.25° increments at all Mach 
numbers of test plan. If radically different spanned lifting surfaces are to be tested later, they should 
also be tested inverted (e.g., variable sweep wings). 


. Level model upright and run baseline configuration at all planned Mach numbers including the same 
-2° to 2° increments. (This need not be a partial run.) 


. Compute model data for @ = 0° (upright run) and @ = 180° (inverted run). 
. Plot Cy vs a for upright configuration and Cy vs -a for inverted run at each Mach number. 
. Prepare upflow table for inclusion in the data reduction program. 


. Recompute upright and inverted data with the upflow correction included to verify that the Cy vs o 
data collapse to a single curve at values of a near 0°. 


. If sideslip testing is included in the test plan, the flow angularity determination is more complicated. 
An empirical interpolation capability has been determined so that flow flow angularity in two planes 
can be estimated at various combinations of model 9 and once 9, and wW, have been determined. 
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DETERMINATION OF UPFLOW ANGLE 


C 
024" 8 
: yf 
\ rs 
eo 17 Model 
\ 9 7 Upright 
N n’ - — — — Inverted 
() 
\ 


2(Quptiow) a * 
a 0.4 


o/ 


UPFLOW ANGLE DETERMINATION 
FOR B SWEEP TESTS AT FIXED a 


1. Make run with model rolled +90°. 
2. Make run with model rolled -90°. 


3. Make analysis similar to that required for a sweep tests. However, 
plot Cy as a function of yaw angle w. (note that w = -B). 


4. 9, is angle correction requide to collapse Cy vs (-B) w data to one 
curve. 
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METHODS FOR DETERMINING SIDEFLOW 


¢ For Aloha Sweep Tests at Fixed Yaw: 
— Determine variation of Cy with yaw (w) (i-e., Cyy) 


— Using data obtained at a = 0° and relevant offset yaw angle solve the 
following equation for w,. 


= Cy (at Wuncorrected* a = 0°) 
ee © ee Wuncorrected* 


(That is, determine adjustment necessary to make Cy = 0 at w= 0°) 


¢ For Beta Sweep Tests at Fixed Alpha: 
— Determine C,,, and Cy, a = 0°. 


— Using data obtained at y = 0° (that is, model rolled +90°) and 
relevant offset yaw angle solve the following equation for W),. 


= e (at Ouncorrecteass B= 0°) - Cys at ) = 
Wy — C Ouncorrected* 
N 


a 


Note: “Uncorrected” here means not corrected for sideflow 
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TO VERIFY SIDEFLOW ANGULARITY 


1. Reprocess data with flow angularity corrections applied. 
2. Plot corrected data to confirm results (See below) 


Model Sweep Fixed 


Corrected Data © Upright a W = 0° 
----- Uncorrected Data © Inverted a wy =0 
@ Upright oO W = constant 
A +90" W a= 0° 
0.8 fe} A IN -90° W a= 0° 
© +90° W a. = constant 
0.12 
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IMPORTANCE OF G@ MEASUREMENTS 


C 
Range = sre — og, 
C JIC, Ww, 


¢ For 
— C-5A, 0.0001 error in C, was equivalent to 1,000 lbs payload 
— HSR vehicle, 0.0001 error in Cp was equivalent to 10,000 Ibs TOGW 


— F-14, error in a of 0.18° resulted in a loss of 200 miles in combat 
radius or 20% of the 500 mile requirement 
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OTHER CORRECTIONS/ADJUSTMENTS 
BASE PRESSURE 
DUCT INTERNAL DRAG 


BASE PRESSURE 


BASE PRESSURE CORRECTIONS 


PURPOSE: To adjust balance measurements for pressure tare force acting on 
model base. This tare force results from model/installation requirements for 
model support and is usually not representative of the pressure conditions 
existing at the base of a real airplane. Therefore, as a matter of convention the 
measured balance forces and moments are routinely adjusted to represent the 
condition of free-stream static pressure existing at the model base. Knowing 
that this adjustment has been made to the model data, it can be adjusted to 
the airplane base configuration if separate tests have been made to determine 
the incremental effects of sting interference and afterbody closure. 


Z 


Model Base 


Model Moment reference 
Moment reference center base center —~r»@ 


Crpase = -IC, pase dA, 
CN base =o J Cp base dA, Pressure force 


Ciitbase = JC, nase LdA 
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BASE PRESSURE CORRECTIONS-Cont. 


Cag - Caieneued 7 Cy base 
Cre = Ciacacuied 7 CN base 
Cas = Criexeued ~ Citase 


Note: Cy, C,, C, can also be corrected for base pressure in a similar way but it is generally not 
necessary since afterbody/installation modifications to the configuration are generally 
symmetric about the vertical plane. 


The integration of pressures is usually done by assigning an area to each base pressure orifice 
and summing the products of the areas and pressures. On some models, there will be orifices in 
the physical base of the model as well as pressure tubes in the open area (often referred to as 
“sting cavity pressures”) at the model base through which the sting passes. (It is not 
recommended that the base or cavity pressures be manifolded to reduce the number of pressure 
transducers required since this can lead to erroneous results. Also, open tubes in the sting cavity 
area are uSually installed with the open area pointed downstream in case there is unexpected 
flow through the model.) 


There is no way to make a base pressure correction if there is flow through the model. 
N 
F, - (Cy pase) AGS re 
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INTERNAL DUCT DRAG 


INTERNAL DRAG CORRECTION 


PURPOSE: To correct balance measurements for tare force due to skin friction 
and momentum losses inside the ducts of flow-through nacelles 


m.Ve + A.(Pe - Po) 


Net Axial Force = [m.V. + A.(Pe - P,)]Cosé,, - mM,V,cosa 


Net Normal Force = m,V,sino - [mM.V. + A.(pe - Po) ]siné, 


Note: ram drag = m,V, 
gross thrust = m,.Ve + Ae(Pe - Do) 


Corrections to the other balance measurements must also be made to account for any cant angle 
of the nacelles. 
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INTERNAL DRAG CORRECTION - Cont. 


Net Drag Force = [m,V, + A,(P. - P,)][Cos®,cosa - sing,sina]— mV, 


For flow-through nacelles 


and 


V> = Mo VyRT, 


Therefore: To determine the internal drag, the following must be measured at the exit of 
the nacelle: 


m, -- Mass flow at exit 
V -- Average flow velocity at exit 
p, -- Average static pressure at exit 
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INTERNAL DRAG CORRECTION - Cont. 


In practice only total pressure profiles and static pressure profiles are actually measured 
since: 


M.VyRT. 
V.=M.VyRT, = — 
6 TEM 2 
Tre = Tro 
M. = f (Pe/Pte) 


Me = PeVeAe = Vy/R 


A.M.V1 + M2 


V Tre 


To make these measurements, a combination total/static pressure rake is used. 


69 


DUCT EXIT PRESSURE SURVEY RAKE NEEDS 


. The rake should be attached to the model to avoid displacement of the rake probes 
relative to the desired duct exit plane locations due to sting/balance deflections during 
testing. 


. Separate tunnel runs are required to obtained internal drag when duct exit plane 
pressure surveys are made. 


. Sufficient number of total pressure probes must be located close to the duct walls to 
adequately define boundary layer. However, care must be taken not to introduce 
blockage at the duct exit plane. 


. Static pressure should be measured by the rake orifices in the duct exit plane or on the 
duct wall at the exit if practical. 


. Care should be taken that the rake support structure be well aft of the duct exit plane 
and not be massive enough to project a pressure field forward to the duct exit plane. A 
duct static pressure tap can give some indication of this occurrence by comparing rake- 
on and rake-off readings. 


. Total and static pressure transducers should be sized and provided with the appropriate 
high quality reference pressure so that good resolution and accuracy is obtained for the 
rake local pressure measurements. 
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OTHER DETERMINATIONS OF INTERNAL DRAG 


¢ Static Calibration Facility (Test facility with flow and force measuring systems) 


Induce duct flow by reducing back pressure (high pressure air supply required) 


Bell-mouth design and fabrication may be required for static test to avoid flow 
separation at inlet lip 


Accurate metering of duct flow required 
Accurate force measuring system required 


Duct wall (near exit) pressures recorded as correlation parameter for use during 
wind tunnel testing 


Boeing has an excellent calibration facility (high costs associated with facility) 


¢ Analytic Computation Methods: 


May be acceptable for some simple duct geometries 
May require duct wall static pressure measurements near duct exit 


Used during HSR testing. However, computational method only computed skin 
friction losses and did not include momentum losses 
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BOUNDARY LAYER TRANSITION 


BOUNDARY LAYER TRANSITION FIXING 


PURPOSE: 


a) Fix transition from laminar to turbulent at known locations on model. 
b) To insure that the boundary layer is essentially turbulent everywhere 
downstream of the trip locations to enable wind tunnel data to be 

extrapolated to flight Reynolds numbers. 


PROCEDURE: 


1. Select location and width of transition trip based on recommendations of 
Peterson in memo to FSRD dated 11/05/65. 

2. Determine critical height of transition trips with recommendations of Braslow 
& Knox in NACA TN 4363, 1958. 

3. For configurations and test conditions where shock induced separated flow 
may occur, consider the recommendations of Blackwell in NASA TN D-5003, 
Jan. 1969. 

4. At low subsonic speeds, there are some who argue that the free transition 
condition in the post stall environment gives more meaningful stability and 
control data. 

5. If using grit, apply to model according to recommendations of Peterson: 

e Narrow roughness band 
e Sparsely distributed grit particles 
e Avoid adhesive buildup during refurbishing. 

6. If using trip dots, follow established procedures. 

7. Check transition strips (grit or trip dots) on a regular basis and refurbish as 
required. 
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Fixed transition 


Wind Tunnel Separation 


Natural transition 
*“— Shock 


“— Separation 
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TRANSITION APPLICATION FOR SUPERCRITICAL FLOWS 


Wind Tunnel - Low R, 


Shock 
Fixed transition forward \ 


Flight - High R, 


Natural transition forward \ / 


Wind Tunnel - Low R, 


Fixed transition rearward 


Reference: NASA TN D-5003 
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MISCELLANEOUS IMPORTANT STUFF 


1. Test Notebooks: 
¢ Good detailed shift notes 
e Leak check records (Be thorough) 
e Describe problems encountered and solutions 
e Sketch devices developed during test (Or you will forget how you did it) 


2. Space In Model: 


e Enough internal room to mount instruments with easy access 
e Large passages to route tubes and leads to instruments or through sting 


3. Support System Space: 

e Large enough openings into sting for tubes, etc. 

e Large enough hole(s) through sting, adapter, knuckles, etc. to route tubes and 
leads 

e Externally routed instrumentation may inhibit ability to roll model or change 
offset hardware 

e Plug sting instrumentation passages with clay of RTV to inhibit flow through 
the sting between the model interior and model support system 
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MISCELLANEOUS IMPORTANT STUFF 2 


4. Special Hardware: 
— Leveling plates should be an integral part of the model design process 


— Consider having the leveling plate weigh the same as any model hardware that is removed to 
level the model (i.e., cover plate) 
— Any special calibration fixtures should also be an integral part of the model design process. This 
would include rigs to calibrate hinge moments of control surfaces. 
5. Zero Returns: 
— Monitor zero returns for instrumentation problems 
* Balance zero shifts 
¢ Plugged pressure orifices 
¢ Slow (crimped or partially plugged) pressures 


6. ESP’s: 
— Large enough routing tubes for Ref., Cal., C1, and C2 pressures 
— Place manifolds for above pressures close to modules 
— Independent check pressures on each ESP (One of two ports each) 
— Monitoring scheme so that ESP recalibration need is obvious on CRT display 


7. Base Pressure Measurements: 


— For performance testing the model base/cavity pressures are as important as force balance data. 
Therefore, it is strongly recommended they be hooked up to individual transducers and not be 
hooked up to ESPs. Therefore, if they are not hooked up to ESPs a decision can be made a any 
time during the test to shut off the ESPs and continue testing. These pressures should be 
measured on very accurate transducers whose power supply can be monitored. 


7 


SOME OFT HEARD STATEMENTS 


¢ We will calibrate after the test is over. 


¢ All we care about are increments. 
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EXHAUST SIMULATION 
TESTING 


Presented by 
Bob Berrier 
4/l17/17 


OUTLINE - PROPULSION TESTING 


¢ Exhaust Simulation Methods 


¢ Model Design Criteria 
— Propulsion Model Decisions 
— Balance Arrangements 
— Support Systems 


¢ Testing Techniques and Data Trends 


EXHAUST SIMULATION METHODS 


EXHAUST SIMULATION METHODS 


TECHNIQUE 


1. Flow Through 


2. Plume Simulator 


ADVANTAGE 


Low cost 
Simulated inlet flow 


Realistic upstream flow field 
for close-coupled propulsion 
systems (transports) 


« Low cost 


Relative low cost 
Safe 


Reasonable simulation of 
initial plume angle 


DISADVANTAGE 


e Flight NPR not matched 
e No thrust measurements 


e Plume blockage effect 
simulated for only one NPR 


e No jet entrainment 
e No thrust measurement 
e Faired-over inlet 


Cold jet - affects downstream 
plume mixing and shape 


Faired-over inlet 


Large volume supply lines 
drives support size up 


No heat transfer tests 


EXHAUST SIMULATION METHODS 


TECHNIQUE 
4. H,O, - 
Hydrogen peroxide 


5. Burners 
Hydrogen 
Methane, etc. 


6. Turbofan Simulators 


7. Real Engine 


ADVANTAGE 


Hot jet - good simulation of cruise 
T,;, plume shape and mixing 


Heat transfer tests 
Small volume supply lines 


Hot jet 
Heat transfer tests 
Small volume supply lines 


Better matching of both inlet and 
exit flow conditions 


Exact match of temperature, 
NPR, plume, chemistry, etc 


DISADVANTAGE 


Cost - Model and H,0, 

Fire hazard, corrosive 
Complex model design 
Balance temperature problems 
Faired-over inlets 


Safety hazard 


Temperature impact on model 
design 
Faired-over inlets 


Does not generally provide 
simultaneous match of inlet and 
exit flows 


Limits model scale 
Complex model, instrumentation 


Requires extensive calibrations in 
special test facility 


Severe limitation on model scale 
Complex 

Cost 

Safety hazard 

Purging of engine exhaust 


RELATION OF ENGINE/NOZZLE/TEST 
PARAMETERS 


ENGINE/NOZZLE OPERATION 


1. Pt,j 
2. Geometry 


C ENGINE AFTERBURNER Test requirements 


Pp. eg] PLA Engine Control] -——=——= } Mach -——— 
ae racemes 


A/B 
Dry : Aircraft air data 
Throttle ( Pilot) sensors 


RPM is a major function of dry PLA; minor function of A/B PLA 
A, is a major function of A/B PLA; minor function of dry PLA 


A, is a major function of A, and Mach; minor function of altitude 


VARIATION OF P, ; 


Max A/B 
Max dry 


Increasing PLA 
—s ¢ Operating NPR schedule 
NPR is different for each 
engine design 
ee Operating NPR is nearly 
a . 
Flow through independent of PLA for 
, max dry and above 
¢ The variation of operating 
RPM NPR with altitude is small 


RPM 
Ay 
Max dry Max A/B 


LA ———=» 


VARIATION OF GEOMETRY 


Increasing 
PLA 


Subsonic Supersonic 


p= 
= 
= 
= 


Max 
A/B 


Part 
A/B 


Max 
Dry 


Part Power 
Cruise 


WW 


Increasing Mach 


NOTES 

1.Nozzle geometry is nearly 
independent of PLA during 
dry power operation 

2.A,increases with PLA 
during A/B operation 

3.A,/ A, increases with Mach 
up to mechanical limits 


4.A,/A, is independent of 
PLA 
but 


since A, varies with PLA 
during A/B operation, A, 
must also vary for A, / A, to 
remain constant 


SUMMARY OF NPR AND GEOMETRY 


By testing Max dry geometry over a range of NPR, the 
entire dry power PLA range can be reasonably 
simulated. 


By testing two throat areas (generally Max dry and Max 
A/B), the entire range of PLA can be covered. 


By testing two expansion ratios for each throat area, 
the entire range of PLA and Mach can be covered. 


A TOTAL OF FOUR NOZZLE CONFIGURATIONS 
TESTED OVER A RANGE OF NPR 


TYPICAL NOZZLE TYPES AND 
PERFORMANCE 


CONVERGENT NOZZLE 


Plume shape with nozzle 


‘ 
sige Underexpanded 
aed ate On design 
ne. || 99 =e= Overexpanded 
Dry 
Desi 
10 NPR. (Pe= Pa) mele 
Design NPR ~ 1.89 
AJA, = 1.0 Choke NPR ~ 1.89 | 
Overexpanded nozzle NPR < Design 
Bee Underexpanded nozzle NPR > Design 
Overexpanded : ; 
nozzle Underexpanded nozzle Most subsonic aircraft 
Pe < Pp e > 00 = 
= have nozzles of this 


1.0 2.0 


NPR type. 


CONVERGENT-DIVERGENT NOZZLE 


Internal separation 


Low s 


Most supersonic aircraft 
have nozzles of this 


type 


Typical internal separation 


SINGLE EXPANSION RAMP NOZZLE 


Notes: 


Upper flap may be used for 
vectoring. 


Lower flap may or may not be 


=== variable. 
“ Free 


Boundary Nozzle may have both an A, ;/A, 
and A.../A, and performance trend 
may have two peaks. 


Since nozzle is not symmetric, 
nozzle normal force and pitching 
will be nonzero, vary with NPR and 


can be quite large. Upper flap 
variation could be used to trim 


these forces and moments. 


Nozzle on B-2 and most hypersonic 
vehicles are of this type. 


1.0 NPR 


MODEL DESIGN CRITERIA 


EXHAUST SIMULATION MODEL 
DESIGN CRITERIA 


Model Design Must Meet Aerodynamic, Propulsive and Structural 
Requirements 


Large Model Size Required to: 


— Accommodate internal propellant lines 


Incorporate flow transfer devices 

Have multiple balances 

Have extensive pressure and temperature measurements 
Properly scaled nozzle details 


Other Model Size Considerations: 
— Subsonic/transonic blockage 
— Shock reflection length considerations 
— Supersonic shock interference 
— Model loads ( Must meet current LHB 1710.3) 


Large Tunnels are Generally Required Because of the 
Size of Powered Models 


PROPULSION AERODYNAMICS RESEARCH 
IS USUALLY A LONG TERM COMMITMENT 


Greater time is involved with complex powered models 
and often requires multiyear funding 


Construction of a complete powered model 
(airframe/balance/support) can take up to one year and 
is very expensive 

Careful assembly, checkout and calibrations are required 
prior to a wind tunnel test and often require a static test 
prior to tunnel entry 

For a six-component balance test, installed balance 
calibrations can last a week or more 

Test programs involve propulsion system variables in 
addition to aerodynamic variables and can easily last 
over one month and require multiple tunnel entries 


PROPULSION MODEL DECISIONS 
¢ Fully or partially metric? 
¢ 0, 1, 2 or more balances? 


¢ Thrust, Drag, and/or Thrust minus Drag 
Measurements? 


TYPICAL MODEL DECISION PROCESS 


Do | need complete aerodynamics ( trim, spillage effects)? 
Do | plan modifications near front of model? 

Do | expect aft-end parametrics to feed well forward? 

Is model close-coupled? 


— Jet/wing interactions, wing-mounted nacelles, 
g g <> lf answer is yes to any 
one of these, then 


jet interactions expected on wing or front of model? 


If answer _is no to all 


FULLY METRIC 


Based on cost, model availability, 


balances, etc 
C4 


PARTIALLY METRIC 


FULL OR PARTIALLY METRIC POWERED MODEL 


Fully Metric Model Partially Metric Model 


Measure complete aerodynamics (trim) Lower balance loads 
All jet interactions included Support interference minimized 


Can make modifications to forward portions of Support interference effects do not have to be 
model determined 


Can be converted to aerodynamic model Fouling minimized 


High balance loads Complete aerodynamics not obtained 


Support interference, particularly for pitch Model modifications usually limited to metric 


May need to determine support interference portion of model 
effects Large internal pressure correction to axial 


Type of metric break seal around strut 
Large internal pressure correction to normal 
and pitch 

Potential for fouling problems 

May cost more 


¢ Thrust - Drag: 


BALANCE ARRANGEMENTS 


SINGLE BALANCE 
¢ Thrust: 
¢ Drag: 


Normally used only for static tests 

Lower balance loads; greater risk for fouling and balance 
dynamics; unrealistic base areas 

Higher balance loads; includes external flow effects on thrust; 
how to separate drag from thrust; lower drag accuracy 


MULTIPLE BALANCES 


¢ Thrust - Drag/Drag: Common 


Thrust/Drag: Common 
Thrust-Drag/Thrust: Not common 
Drag/Thrust - Drag: Not common 


Drag/Drag: 


Not common 


PROPULSION MODEL SUPPORT 
SYSTEMS 


PROPULSION MODEL SUPPORTS 


1.Same as aero models 
a)Loads 
b) Divergence 
c) Test section position 
d)Alpha/beta requirements 


e) Instrumentation leads MINIMUM 
= 
2. Must allow geometry SUPPORT 
simulation in region of 
interest (inlet/nozzle) INTERFERENCE 


+ 


3. Must supply exhaust 
simulation to model 


FOR PROPULSION MODELS, SUPPORT SYSTEMS 
ARE A COMPROMISE BETWEEN REQUIREMENTS’ 


TECHNIQUE 
1. Sting 


2. Annular-flow sting 


3. Solid plume 
simulators 


4. Strut/sting-strut 


5. Bifurcated 
(wing-tip support) 


ADVANTAGE 


e Minimum support interference 


No aft-end distortion, particularly 
for A/B 


Small support interference 


Minimum support interference 
Simple/low cost 


No aft-end distortion 


Exhaust medium supplied to front 
end of model 


No aft-end distortion 


Exhaust medium supplied to front 
end of model 


Can be used to simulate wings 


MODEL SUPPORT TYPES 


DISADVANTAGE 


e Distorts aft-end geometry 


Complex arrangement when 
thrust is measured 


Distorts plume shape 

Large sting to dry throat ratio 
Plume blockage effect simulated 
at only one NPR 

No jet entrainment 

No thrust measurements 


Support interference on aft-end 


Significant interference on 
complete metric model 


Struts (wings) distorted at tip 
Doubles isolated model 
interference 

Complete metric model very 
difficult 


MODEL SUPPORT TYPES 


TECHNIQUE 
6. Extended forebody 


7. Semispan from 
wall 


8. Semispan model 


ADVANTAGE 


No aft-end distortion 


Exhaust medium supplied to front 
end of model 


Minimum direct support 
interference 


Larger model scale possible 
No aft-end distortion 


Ease of bridging balance with 
exhaust and instrumentation lines 


No aft-end distortion 
Ease of bridging balance 
Support simulates wing 
Real forebody effects 


DISADVANTAGE 


Unrealistic boundary layer 
Additional complexity to simulate 
boundary layer 


Distorted forebody for complete 
metric models 


Wall interference and boundary 
layer distortions 

No yaw capability 

Special balances with large yaw 
component required 
Nocross-talk about centerline 


Unsymmetric wing distortion 
Special balance 

Axial metric break and seal 
Unknowns ( not been done) 


MODEL SUPPORT VARIATIONS 


FULLY METRIC MODELS PARTIALLY METRIC MODELS 


Forward Strut - Single Balance 


Wing-Tip Suppor 


Forward Strut - Two Balance 


STING-SUPPORT INSTALLATIONS 


STING-STRUT INSTALLATIONS 


WING-TIP SUPPORT INSTALLATIONS 


VARIOUS OTHER INSTALLATIONS 


DETERMINE STING SUPPORT INTERFERENCE 


¢ Normal Run 
— Model force or moment 


— Interference of sting on bottom of 


Live Sting 
model 


Top Strut Live 


AES ESE ERE BREE RRR ETI - Tare RunA 
(ES a — Interference of top strut on model 


OI = _ interference Gt sting on bottom at 
Dummy Stig adel 


¢ Tare Run B 
a — Interference of top strut on model 


Interference of sting = Tare Run A - Tare Run B 


Top Strut Live 


Corrected Force = Normal Run - (Tare Run A - Tare Run B) 


DETERMINE STRUT SUPPORT INTERFERENCE 


¢ Normal Run 
— Model force 


— Interference of bottom strut on 
model 


Bottom Strut Live 


Top Strut Live 


¢ Tare RunA 
—Interference of top strut on 
model 
—Interference of bottom strut on 
model 


Bottom Strut Dummy 


Top Strut Live 


e Tare Run B 


— Interference of top strut on 
model 


Interference of bottom strut = Tare Run A- Tare Run B 
Corrected Force = Normal Run - (Tare Run A - Tare Run B) 
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DATA REDUCTION REQUIREMENTS 
FOR PROPULSION TESTING 


Presented by 
Larry Leavitt 
4/17/17 


12/13/06 


OUTLINE 


* Propulsion Tares 
— Propulsion Simulation Systems 


— From Whence We Came, The ‘Good Old Days’ 
— Propulsion Simulation Systems Using Twin Bellows 
— Instrumentation/Hardware Requirements 
= Accurate Flow Measurement Apparatus 
= Calibration Hardware 
— Where Do Tares Come From? 
— How Do We Find Them? 
— Description of Various Tares 
= High Restraints 
= Jet-Off Force/Moment Interactions 
= Loading/Blowing Interactions 
# Axial Momentum Tares 
¢ Other Data Reduction Requirements 
— Nozzle Flow Parameters 
— Pressure Coefficients & Integrated Forces 
— Aerodynamic Analysis Techniques (Thrust Removal) 
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PROPULSION TARES 


PROPULSION SIMULATION SYSTEMS 


Twin Engine/Metal Bellows 


Twin Engine With Metal Bellows 
And Plastic Diaphragms 


. 
Single Engine/Metal Bellows 


System Using Variation of Long Pipe 


Secondary plenum 


Structural support cart 
(removed for clarity) —_ | 


Secondary S-tubes ——— 


Primary S-tubes — 


Dual Flow Test Stand/S-Tubes 
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Things That Haven’t Worked Well for Us and Other 
Considerations 


¢ Aijr transfer systems that bridge metric/non-metric break using hardware that 
slide or rotate on an o-ring(s) (hysteresis a problem) 

¢ Piston type force balance with labyrinth seals — you modified pressure to 
remove a delta pressure across the system. Knowing area, one could 
calculate the force. 

¢ Relying on loose high pressure hoses to bridge the metric/non-metric break 
(hysteresis and other uncertainties like strut movement) 

¢ We never had much luck with flow thru force balances (Others have. 16TT 
temperatures may have been the problem) 

¢ Some systems are very susceptible to differences in external flow 
temperature and internal flows (propulsion simulation air). 


Note: any system that produces gradients of different temperatures across 
different portions of the force balance give problems. 

Goal: A repeatable (yet flexible) means of getting high pressure air (or other 
propulsion simulation medium) across the metric/non-metric break 


FROM WHENCE WE CAME! 


° eally Days of Propulsion Testing ( Late 1960s) 


Original air-powered system mounted on a fixed strut (from tunnel floor) was a replacement 
for a H,O, hot-jet system that was put into operation in early 1950s 


— Could only test at zero angle of attack 
— Only momentum corrections were required 


— Weight flow measured with rotary type flow meters (low accuracy because frequency 
converter needed to determine weight flow) 


¢ Mid Technology (Mid 1970s) 
Model mounted on a sting-strut that provided angle of attack variation 
— Realization that blowing/loading tares important with the advent of vectored thrust research 
— Found during testing of new twin-jet powered model designed for thrust vectoring 
— Mostly 3 component testing 


— Multiple Critical Flow Venturi System put in operation with the capability to measure flows 
from 1 |b/sec to 40 Ibs/sec. This system depends only on a single pressure measurement. 


¢ Latest and Greatest (1980s - ?) 


— Continuous improvements and refinements as testing requirements increased 
— 6-component corrections necessary for testing multi-axis thrust vectoring nozzles 
— Automated data reduction technique for determining tare constants put in operation 
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PROPULSION SIMULATION SYSTEMS WITH 
BELLOWS FLOW TRANSFER DEVICES 
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SINGLE-ENGINE PROPULSION SIMULATION SYSTEM 


—s}$ 


<E 


Flow transfer device 
(bellows system) 


¢ Workhorse propulsion simulation system developed in late 1960s to replace 
HO, hot jet system 
¢ Featured a unique, compact, twin-bellows flow transfer device designed to 


eliminate momentum tare forces 


¢ Used continuously in the 16-Foot Transonic Tunnel and Jet-Exit Test Facility 
from 1968 to 2005. 
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TWIN-ENGINE PROPULSION SIMULATION SYSTEM 


Clearance holes for sonic nozzles 


Sonic nozzles 
Center ring 


Metal bellows 


Airflow 
SEE el 


PCH (psi) 
Flange 


Cavity vented to model 


PBEL/(psi) sleeve internal pressure 
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INSTRUMENTATION/HARDWARE 
REQUIREMENTS 


A VERY ACCURATE FLOW MEASURING SYSTEM 
REQUIRED 


Multiple Critical Flow Venturi System 
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STRATFORD CALIBRATION NOZZLES 
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CALIBRATION HARDWARE NEEDS TO BE PART OF BASIC 
MODEL DESIGN 


Twin Jet Model with Stratford Nozzles 


& = Sa 


F-18 Model with Calibration Fixture 13 
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WHERE DO TARES COME FROM? 
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TWO TYPE TARES NEED TO BE DETERMINED 


¢ JET-OFF TARES 


— Changes in the spring constant of the balance resulting from the Flow Transfer Device 
(bellows assembly), steel tubing, cables, etc. bridging the metric and nonmetric 
portions of the model 


Ketlows 


Balance 


K 


tubes, cables, etc 


¢ JET-ON TARES 
— Changes in the spring constant of the bellows resulting from pressurization of the 
bellows, local pressure differences from one end of the bellows to the other, and 
misalignment of the momentum transfer vector 


A 


pressure/momentum 


Ktubes, cables, etc 15 
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MOMENTUM TARE FORCES 


Pressure differential between 
bellows ends 


Supply air 


Misaligned vector 
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HOW DO YOU FIND THEM? 


THREE THINGS YOU MUST KNOW OR DETERMINE 


e Tare Free Balance Characteristics 
— Bare balance calibration 


¢ Effect of Restraints Across Balance at Jet-off Conditions: 


— Loading calibration of balance with fully assembled model 
including all tubing, cables, etc. 


Effect of Bellows Pressurization and ‘Flowing’ Restraints 
Across Balance at Jet-on Conditions: 


— Loading/blowing calibration of balance with fully assembled 
model including all tubing, cables, etc. 
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PROPULSION TARES CALIBRATION PROCEDURE 


¢ System Leak Check -- NO Leak Is a Good Leak! 


¢ Perform Jet-Off Balance Loadings 
— All components of interest 


¢ Perform Jet-On Balance Loadings 
— All components of interest 
— Over expected range of loads, bellows pressure and throat areas 


¢ Perform Axial Momentum Tare Runs 
— Over expected range of throat areas 


F. J. Capone - 12/13/06 v 


TYPICAL BALANCE LOADINGS CONSIDERATIONS 


¢ Provide as many loads possible by hanging because pulleys often 
provide hysterisis problems 


¢ Loadings should be carefully directed through the balance moment 
reference center. Be sure your loading fixtures are properly placed 
and model is level at each point (unless pulleys are attached to the 
model support system) 


¢ Hang all loads to expected maximum values. If expected maximum 
values are unknown, load to the balance limits. (Positive and 
negative loadings were suggested for the automated tare program) 
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TYPICAL LOADING/BLOWING MATRIX 


Example shown for Normal loadings 


_Normal_| __NPR___| Cal Nozzle Area 
[0 | 1.0 (JET-OFF) 
[200 fp 


¢ As needed load: 
+ Pitch, = Roll, == Yaw, = Side 


¢ Typically pitch, roll and yaw are 
loaded by transfer of weights. For 


example, would have constant NF 
while PM varies 

¢ Roll tares may be questionable 
because bellows do like to be 
twisted 
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DESCRIPTION OF VARIOUS TARES 
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HIGH RESTRAINTS 


Measured 


Kretiows 


Balance 


K 


tubes, cables, etc 


Measured sensitivity constant 


HIR, = Balance sensitivity constant ; 


High restraint constants account for the loss in component sensitivity 
due to additional “spring constants” across the balance. 


For a given load, millivolt output will be less for the installed balance 


than for the bare balance. 
High restraints constants SHOULD ALWAYS BE POSITIVE 
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JET-OFF FORCE/MOMENT INTERACTIONS 


Example: Determine Normal Force Tares for 6 Component Test 


fap) (99) 
So) = al 9 oo) 
: 3 & = b 
a - a a 
NFicad NF cad NF iad NF iad NF iad 


Where: AAF1,3 =AF1,3xricag - AF 1,3ye-o, APM1,3 = PM1,3ycicag - PM1,3yp=o, ete. 


Compute jet-off runs with high restraint constants and XK’s = 0 and 
KMOM = 2 to determine tares. 


For a 3 component test, need to look at effects of: 
— NF onAF and PM 
— PMonAF and NF 
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LOADING/BLOWING INTERACTIONS 


Example: Determine Normal Force/Air Effects on NF Tare for 
Particular Calibration Nozzle 


NF Load 
800 
ee ANF 1,3 7,1 
se) J = xk4 ? 
L oe APBEL ae 
z 200 
0 xk46,1 for particular 
static run only 
APBEL NF Load 


where ANF1,3 = NF1,3yer on - NF1,3yer ogr » ete. 


Compute all runs with high restraints and jet-off constants with correct weight flow. 
The effects of NF/Air on DPM, DRM, DYM, DSF must also be determined. 


Likewise, the effects of PM/Air, RM/Air, YM/Air, and SF/Air on all balance must be 
computed as required 


Recompute data, CHECK RESULTS, iterate 
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Typical Tare Equation 


Jet off tares 


TAREN = XK(1,1) + XK(2,1)*FN + XK(3,1)*PM+ XK(4,1)*RM + 
XK(5,1)*YM + XK(6,1)*SF + DELP {(XK(46,1) + XK(47,1)*FN + 
XK(48,1)*PM + XK(49,1)*RM + XK(50,1)*YM + XK(51,1)*SF + 
FNO*XK(52,1) + FNO2*XK(53,1) + PMO*XK(54,1) + PMO2*XK(55,1) + 
AREA*[XK(56,1) + XK(57,1)*EN + XK(58,1)*PM + XK(59,1)*RM + 
XK(60,1)*YM + XK(61,1)*SF] + AREA2 [XK(62,1) + XK(63,1)*FN + 
XK(64,1)*PM + XK(65,1)*RM + XK(66,1)*YM + XK(67,1)*SF]} 


i Blowing tares 


DELP — Bellows Pressure minus free-stream static pressure 
FNO — Normal force preload 

PMO -— Pitching moment preload 

AREA — Nozzle area 
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AXIAL MOMENTUM TARES 


¢ Determined from the zero-load calibration runs 

¢ All high restraint and all nonblowing corrections applied to data 

¢ Mass flow must be correct, i.e., nozzle discharge coefficient is correct 
¢ Determined by: FAMOM = AF(1,3) - (FJCON/FI)*FI 


¢ FJCON/FI is found by table look-up from historical data on Stratford choke 
nozzles. Any convergent nozzle could be used as a calibration nozzle as 
long as its performance is known. 


DPEB 
0 80 120 


Cal Nozzle 2 


As can be seen, FAMOM is a 
function of both bellows pressure 
and nozzle throat area 


Cal Nozzle 3 
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OTHER DATA REDUCTION REQUIREMENTS 
FOR PROPULSION TESTING 
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INTEGRATED DATA REDUCTION PROGRAM 


Tunnel Parameters 


Balance Loads & 


Nozzle Flow 
Parameters Model Attitude 

Skin Friction Pressure Coefficients & 
Calculations Integrated Forces 


Aerodynamic Analyses 
(Thrust Removal) 
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NOZZLE FLOW PARAMETERS 


Nozzle Exhaust Conditions 


— Measurement of average total 
pressure and temperature usually 
made at a charging station located 
upstream of nozzle exit 


Nozzle Pressure Ratio 


— Ratio of average jet total pressure to 
free-stream static pressure 


— NPR = p,;/p,, 


Nozzle Mass Flow Se es oe 
— Total mass flow w, obtained from NOZzlé€ ideal | Nrust 


venturi measurements — Ideal Thrust = w; = f(NPR, Tj, Wp) 
— Ideal mass flow w, = f(p;;, Ty)) 


RT., Y. 
+1 F, = w, a 2y- 1- f4 
¢ \y-1}| \NPR 


— Nozzle performance parameter, F/F; 


Pr gy ( 2 


Fy dR 


— Discharge coefficient = w, / w; 
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PRESSURE COEFFICIENTS & INTEGRATED FORCES 


¢ Eight groups of pressure coefficients could be computed in this module with 
a maximum of 125 measurements per group 
¢ Calculated quantities from each group consisted of: 
— C, = (p - p.,)/., 
— CE = D(C, Ay ! Avert 
~ Crip = X(C, A; ! Aver 
= Crip = X(C, A ! Aver Cref 
— Cyp = X(Cp A / Aver Cret 
~~ Crp = X(C, A; ! Aves Cref 
= Cy, = D(C, Ay ! Avet 
— Cam = U(CpiA/ Aver Cret 
— Cp, = C-,cosa + Cy sina 
— C_,=Cy,pcosa Cr,sina 
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AERODYNAMIC ANALYSES TECHNIQUES 
(THRUST REMOVAL) 


Several options were available to remove thrust and to obtain the various 
aerodynamic and aeropropulsive parameters required for data analyses. These 
options worked for both fully and partially metric models. Computed inputs from 
the previously shown modules were required. 


¢ Some Simplified Thrust Removal Equations 


Computed jet axial force -- C-j. = p,/q., [k(NPR) + I] 

Thrust removed axial force coefficient -- Ca sero = Ca - Cesc 
Computed jet normal force -- Cyj¢ = p,/q., [kK(NPR) + 1] 
Thrust removed normal force coefficient -- Cy sero = Cn - Croc 
Thrust vector angle -- 6, = Cyyo / Ceyc 

Constants k and i are determined from static runs 

Similar calculations can be made for the forces and moments 
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